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Abstract 

Terahertz  radiation  is  an  emerging  field  that  has  far  reaching 
applications.  There  is  a  need  for  portable  and  affordable  terahertz  sources  that 
may  be  used  for  detection  of  structural  weaknesses  in  aerospace  composites. 

One  possibility  for  terahertz  generation  is  a  femtosecond  Cr:LiSAF  laser. 
However,  controlling  intracavity  dispersion  is  necessary  in  order  to  get  the 
ultrashort  pulses  needed  to  generate  terahertz  radiation.  This  research  measures 
dispersion  of  the  chirped  femtosecond  intracavity  mirrors  by  employing  white- 
light  interferometry  and  explains  the  theory  behind  the  curve-fitting  process 
used  to  calculate  dispersion.  To  compensate  for  the  narrow  reflectivity  band  of 
the  mirrors,  a  CaF2  (calcium  fluoride)  window  was  used  to  increase  the  number 
of  spectral  fringes  in  the  interferogram.  A  cavity  was  designed  using  these 
mirrors  and  the  Cr:LiSAF  laser  was  successfully  modelocked. 
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MEASURING  THE  DISPERSION  OF  LASER  CAVITY  MIRRORS  USING 


WHITE-LIGHT  INTERFEROMETRY 

I.  Introduction 

Modelocked,  diode-pumped  CrrLiSAF  (chromium-doped  lithium  strontium 
aluminum  fluoride)  lasers  offer  a  portable  and  affordable  alternative  to 
commercial  Ti:Sapphire  (titanium  sapphire)  laser  systems.  When  modelocked, 
they  are  energy  efficient  and  capable  of  producing  femtosecond  pulses  (1015  s  = 

1  fs)  that  can  be  used  to  generate  terahertz  radiation.  In  addition,  Cr:LiSAF 
lasers  can  be  used  in  detection  devices  to  locate  weaknesses  in  aerospace 
composites  and  to  observe  ultrafast  processes  such  as  a  molecule  undergoing  a 
chemical  reaction. 

Controlling  cavity  dispersion  is  crucial  to  modelocking  a  laser.  This 
research  used  white-light  interferometry  to  measure  dispersion  in  laser  cavity 
mirrors.  The  dispersion  results  were  then  used  to  engineer  a  modelocked  laser 
cavity  with  the  goal  of  generating  terahertz  radiation. 

Dispersion  occurs  when  the  phase  velocity  of  a  wave  depends  on  its 
frequency.  For  optical  fields,  the  angular  frequency  is  u  =  c/ k() .  where  c  is  the 
speed  of  light  in  a  vacuum  and  kQ  =  27t/Aq  .  For  a  pulse,  the  phase  velocity  is 
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vp  =  w/k  ,  the  group  velocity  is  1 Jv  =  dk/dcu ,  and  the  group-velocity  dispersion 
is  GVD  =  d2k/ div2 .  The  main  concern  of  this  paper  is  group-delay  dispersion 
GDD = GVD  x  L  ,  which  accounts  for  the  length  of  dispersive  media,  L,  and  also 
applies  to  mirrors  and  laser  cavities.  It  is  commonly  described  in  units  of  fs2. 
Positive  GDD  causes  a  pulse  to  broaden  temporally,  making  it  difficult  to 
achieve  the  desired  effects  of  ultrashort  pulses.  This  research  focused  on  optical 
elements  that  contribute  negative  GDD  and  counteract  pulse  spreading. 

This  experiment  is  a  continuation  of  the  thesis  work  done  by  David 
Jones,  a  previous  AFIT  student  [1],  The  group  delay  dispersion  of  the  optical 
components  in  the  Cr:LiSAF  laser  cavity  were  measured  using  white-light 
interferometry.  This  information  was  used  to  analyze  dispersion  in  the  laser 
cavity.  The  goal  was  to  create  negative  GDD,  a  condition  for  solit-on-like 
modelocking.  This  concept  is  expanded  in  Appendix  A.  Derivation  of  the  Non- 
Linear  Schrodinger  Equation.  Characterization  of  the  CW  beam  and 
modelocking  the  laser  precede  the  goal  of  generating  terahertz  radiation. 
Background 

Reliable  sources  and  detectors  of  terahertz  radiation  (T-rays)  were  hard 
to  find  until  the  1990s.  Now  they  are  available  on  a  limited  basis,  typically  large 
systems  for  scientific  research.  In  2006,  electrical  engineer  Daniel  Mittleman, 
from  Rice  University's  T-ray  lab  said:  “The  technology  is  very  young. 

Terahertz  is  now  where  X-rays  were  in  1905  [2].” 
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T-rays  can  be  used  for  detection  in  organic  and  inorganic  materials, 


communication,  and  astronomy.  They  are  capable  of  penetrating  non¬ 
conducting  materials  like  clothing  and  plastic,  but  not  metal  or  water.  Unlike 
X-rays,  T-rays  do  not  damage  body  tissue,  so  they  are  safe  for  medical  imaging. 
When  looking  at  tissue,  T-rays  are  able  to  detect  differences  in  density  that 
make  them  useful  for  finding  tumors. 

T-rays  have  astronomical  applications  as  well.  The  Herschel  Space 
Observatory,  a  satellite  due  to  launch  in  2008,  is  the  terahertz  version  of  the 
Hubble  telescope  [2],  DOD  interest  includes  non-destructive  testing  of  aerospace 
composites,  T-ray  radar,  and  secure  communications.  Since  many  aircraft 
structures  are  built  from  composite  materials,  it  is  important  to  identify  small 
cracks  and  de laminations.  With  a  small,  portable,  T-ray  source,  this  and  many 
other  applications  would  be  possible. 

Previous  Work 

All-solid-state  laser  systems  became  available  over  a  decade  ago.  Since 
then,  the  use  of  tunable  femtosecond  lasers  has  grown  to  the  point  where  they 
have  become  a  fundamental  tool  in  many  areas  of  science,  including  physics, 
chemistry,  and  biology  [3].  Materials  scientists  probe  the  dynamics  of  electrons 
in  carbon  nanotubes  by  carefully  tuning  a  laser  pulse’s  wavelength  for  time- 
resolved  ultrafast  spectroscopy.  Physicists  also  use  tunable  ultrafast  laser 
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sources  to  explore  the  specific  responses  of  different  nonlinear  optical  materials 
[3], 

Cr:LiSAF  lasers  were  first  demonstrated  in  the  1990s  by  M.  J.  P.  Dymott 
and  A.  I.  Ferguson  [4],  Their  goal  was  to  generate  short  pulses  with  a  high 
repetition  rate  and  high  power.  Their  motivation  was  to  create  an  all-solid-state 
alternative  to  the  large  TirSapphire  laser  [5].  All-solid-state  lasers  are  logistically 
simpler  than  gas  lasers  or  dye  lasers;  and  the  Cr:LiSAF  laser  is  a  nice  alternative 
to  the  Ti:Sapphire  laser  because  it  has  a  smaller  footprint  and  an  uncomplicated 
laser  cavity. 

A  recently  published  experiment  used  spectral  filters  on  flash  lamp  pumps 
to  reduce  waste  heat  in  the  laser  crystal.  The  Centro  de  Lasers  e  Apligacoes  in 
Sao  Paulo,  Brazil  reported  capabilities  of  30  W  average  power  at  851  nm,  with  2 
J  pulses  at  15  Hz.  With  a  more  restrictive  filter,  they  showed  average  powers  of 
20  W  but  peak  powers  of  lOkW  at  30  Hz  [6].  This  group  focused  on  increasing 
power  from  Cr:LiSAF  lasers,  not  T-ray  production. 

Increasing  continuous  wave  (CW)  power  is  worthwhile,  but  short  pulses 
inherently  have  high  peak  powers,  so  more  can  be  done  with  a  lower  average 
power.  The  intensity  of  a  CW  laser,  I  ,  is  simply  defined  as 

_  Power 
cw  Area  ' 
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The  intensity  of  a  pulse,  however,  must  consider  the  pulsewidth,  Tp ,  and  the 
time  between  pulses,  trt  ,  which  is  derived  from  the  cavity  length,  L,  using 


1  2  L 


T  RT 

V  C 


(1.2) 


where  v  is  the  pulse  frequency,  and  c  is  the  speed  of  light.  Taking  the  unique 
qualities  of  a  pulse  into  consideration  yields: 


P  X  T 

p  _  CW  RT 


(1.3) 


V 

Compared  to  trt  ,  r  is  very  short,  so  one  can  see  in  (1.3)  that  peak  power  can 
be  significant,  even  if  the  averaged  CW  power  is  unremarkable  [7].  Modelocking 
is  the  easiest  way  to  get  short  pulses  from  a  laser. 

Modelocking  can  only  be  achieved  if  cavity  dispersion  is  controlled. 
Previous  work  on  measuring  dispersive  media,  both  reflective  and  transmissive, 
used  white-light  interferometry.  That  experimental  technique  is  very  different 
than  the  work  presented  in  this  research.  It  used  a  compensated  Michelson 
interferometer  that  is  later  unbalanced  with  the  addition  of  dispersive  media.  A 
He-Ne  laser  is  co-aligned  with  the  collimated  white-light  beam,  and  the  fringes 
are  sampled  using  a  second  detector.  By  linking  this  detector  to  a  piezoelectric 
translating  stage,  the  movement  of  the  delay  arm  is  measured  to  sub-wavelength 
accuracy.  To  find  dispersion,  a  reference  is  first  taken  without  the  dispersive 
media  using  a  square-law  detector  that  sees  the  intensity  of  the  spatial  fringes  as 
the  delay  arm  is  scanned.  A  second  collection  is  taken  with  the  dispersive 
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media,  and  then  normalized  using  the  reference  data.  The  resulting  plot- 
contains  the  phase  information.  By  fitting  a  fifth-degree  polynomial  to  the 
curve,  it  is  possible  to  discover  the  dispersive  properties  of  the  media  with  less 
than  4%  error  [8-10]. 

Research  Objective 

Modelocked  lasing  from  this  cavity  depends  on  compensating  for  the 
positive  GDD  introduced  by  the  Cr:LiSAF  crystal.  White-light  interferometry 
was  used  to  measure  the  amount  of  GDD  contributed  by  mirrors  in  the  cavity. 
Chapter  2  covers  GDD  in  further  detail. 

This  research  furthered  the  development  of  a  terahertz  radiation  source 
that  is  small,  inexpensive,  portable,  and  rugged.  The  concepts  used  in  this 
system  may  be  miniaturized,  or  expanded  into  a  larger  but  still  simple  system. 
The  applications  of  this  technology  are  far  reaching  -  including  biomedical 
imaging,  aerospace  composite  testing,  airport  security  systems,  and 
communications. 
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II.  Theory 


This  section  covers  the  Cr:LiSAF  laser  cavity  design,  including  selection 
of  the  gain  medium,  necessary  intracavity  corrections,  dispersion  control,  and 
modelocking.  It  also  covers  the  mathematical  theory  supporting  the  use  of  a 
Michelson  interferometer  to  measure  GDD  in  the  laser  cavity  mirrors.  The  MI 
theory  is  extended  to  explain  the  curve-fitting  procedure  used  for  data  analysis. 
Finally,  Grenouille  FROG,  another  method  of  GDD  measurement  is  described. 
Cavity  Design 

The  basic  cavity  design  (Figure  1)  has  the  following  components:  a 
polarizing  beam  splitter  (PBS);  two  spherical  mirrors  (SM)  (R  =  99.9%);  two 
planar  mirrors;  one  with  high  reflectivity  (99.9%);  and  one  with  lower  reflectivity 
(99.0%)  that  acts  as  an  output  coupler  (OC).  A  3  mm  long  x  3  mm  diameter 
Cr:LiSAF  crystal  with  5%  doping  lies  between  the  two  spherical  mirrors.  Two 
AlGalnP  (aluminum  gallium  indium  phosphide)  diode  lasers  are  aligned  such 
that  one  is  polarized  vertically  while  one  is  polarized  horizontally  and  are 
combined  in  a  PBS.  The  combined  beam  passes  through  a  lens  (f  =  5  cm), 
through  the  back  of  a  spherical  mirror,  and  focuses  on  the  crystal.  The  lens 
focuses  the  beam  so  the  volume  of  the  cavity  mode  matches  the  gain  volume  of 
the  crystal.  To  form  a  ring  cavity  the  two  planar  mirrors  are  adjusted  to  reflect 
the  beam  from  each  spherical  mirror  onto  the  opposite  planar  mirror.  The 
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spherical  mirrors  are  adjusted  so  that-  the  beam  waist  is  halfway  between  the 


planar  mirrors.  In  this  configuration  the  CW  beam  propagates  bi-directionally. 
However,  when  the  laser  is  modelocked,  unidirectional  propagation  is  preferred. 
This  effect  is  due  to  gain  competition  and  the  nonlinear  loss  mechanism  induced 
via  the  Kerr  lens.  When  a  saturable  absorber  mirror  (SAM)  is  introduced  into 
the  cavity,  the  beam  path  is  linear  (Figure  2). 

The  Cr:LiSAF  crystal  was  selected  because  of  its  absorption  and  emission 
spectra  (Figure  3).  There  is  a  broad  absorption  band  centered  near  650  nm  and 
an  emission  spectrum  centered  near  840  nm.  Semiconductor  laser  diodes,  made 
with  AlGalnP,  operate  at  this  peak  absorption  wavelength  and  are  readily 
available  at  a  relatively  low  cost.  These  diode  lasers  are  the  same  lasers  used  in 
DVD  burners.  Their  beam  quality  is  excellent  since  burners  require  the  beam  to 
focus  to  a  small  spot  size. 


8 


PBS  Lens  SM  Cr:LiSAF  SM 


Figure  1.  Cavity  design  and  laser  path  for  linear  (red)  and  ring  (red  &  blue) 


PBS  Lens  CVI-H  Cr:LiSAF  CVI-H 
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The  emission  spectrum  of  Cr:LiSAF  is  similar  to  the  more  common 


commercial  Ti:Sapphire  lasers,  making  it  easy  to  find  chirped  mirrors,  and  other 
optical  components.  Cr:LiSAF  is  an  excellent  gain  material  with  good  energy 
storage  and  slope  efficiency,  which  makes  it  ideal  for  use  in  pulsed  lasers  with 
high  peak  powers.  So,  using  this  crystal  as  the  gain  medium  allows  for  an 
efficient,  cost-effective  device. 

6 

4 
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Cavity  Correction 

There  are  two  types  of  correction  needed  in  the  cavity:  astigmatism  and 
GVD.  Astigmatism  is  introduced  into  the  beam  as  it  passes  through  the 
Brewster’s  angle  cut  laser  crystal.  The  vertical  and  horizontal  planes  of  the 
beam  follow  different  paths,  resulting  in  an  elliptically  shaped  beam  hitting  the 
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Figure  3.  Cr:LiSAF  absorption  and  emission  spectra  [11]. 
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spherical  mirrors.  With  a  small  angle  of  incidence,  the  spherical  mirrors  correct 


this  astigmatism,  restoring  the  beam  to  a  circular  cross-section. 

GVD  refers  to  a  specific  effect  of  an  optical  medium  on  a  pulse  of  light. 
Pulse  envelopes  are  made  up  of  a  number  of  frequency  components,  which  are 
added  together  coherently.  In  a  dispersive  system,  these  frequency  components 
travel  at  different  velocities,  so  the  initially  short  pulse  may  broaden  [12]. 
Positive  GVD  slows  the  shorter  wavelengths  with  respect  to  the  longer 
wavelengths  (Figure  4).  This  delay  results  in  temporally  broadened  pulses. 


Delayed  &  broadened 
pulse 


R 


B 


Figure  4.  Positive  GVD  for  a  pulse  in  a  dispersive  medium  [13]. 


The  Sellmeier  equation  can  be  used  to  find  the  amount  of  dispersion  in 
terms  of  wavelength  for  a  given  material.  From  Diels,  GVD  is  250  fs2/cm  for  2% 
doped  Cr:LiSAF.  In  this  3  mm  long  5%  doped  crystal,  the  GDD  is  75  fs2  [10]. 

In  a  linear  cavity,  the  beam  passes  through  the  crystal  twice,  so  the  GDD 
doubles. 

Chirped  mirrors  have  special  multilayer  dielectric  coatings  with  periodic 
spatial  frequency.  They  reflect  shorter  wavelengths  at  the  layers  closest  to  the 
front  surface,  while  allowing  longer  wavelengths  to  penetrate  deeper  layers  before 
being  reflected  (Figure  5). 
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Figure  5.  Operation  principle  of  a  chirped  mirror.  Long  wavelengths  penetrate 
deeper  into  the  mirror  structure  and  thus  experience  a  delay  [14]. 

This  causes  the  delay  of  the  longer  wavelengths  needed  to  create  negative 
GDD.  It  counters  the  effect  of  the  crystal  and  shortens  the  pulse.  Chirped 
mirrors  are  mainly  used  inside  ThSapphire  oscillators  that  produce  pulses  shorter 
than  30  fs  [3].  The  chirped  mirrors  in  this  experiment  are  optimized  for 
femtosecond-pulsed  laser  systems.  The  advantage  of  chirped  mirrors  over  other 
types  of  dispersion  compensation  is  negative  GDD  over  the  whole  bandwidth  of 
femtosecond  pulses.  Assuming  a  Gaussian  pulse,  the  negative  GDD  results  in  an 
ideal  linear  chirp  [15]. 

Chirp  is  the  name  for  a  change  in  frequency  within  a  pulse  envelope.  As 
a  pulse  propagates  through  the  crystal,  it  undergoes  self-phase  modulation 
(SPM).  This  is  a  small  non-linear  effect  where  the  pulse  experiences  a  phase 
shift,  and  subsequent  change  in  the  frequency  spectrum  (Figure  6).  The 
frequency  either  increases  (up  chirped)  or  decreases  (down  chirped)  with  time. 

An  up  chirped  pulse  is  shown  in  Figure  7.  SPM  broadens  the  spectral 
bandwidth,  and  gives  short  pulses  in  the  time  domain  via  the  Fourier  transform, 
as  there  is  an  inverse  relationship  between  bandwidth  and  temporal  pulse  length. 
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Figure  6.  A  pulse  (top  curve)  propagating  through  a  non-linear  medium 
undergoes  a  self-frequency  shift  (bottom  curve)  due  to  SPM.  The  front  of  the 
pulse  is  shifted  to  lower  frequencies,  the  back  to  higher  frequencies.  In  the 
center  of  the  pulse,  the  frequency  shift  is  approximately  linear  [10]. 


Figure  7.  A  positively  chirped  pulse  in  the  time  domain,  resulting  from  SPM. 
In  this  figure,  the  amplitude  is  a  Gaussian  function.  The  phase  is  quadratic, 
resulting  in  an  instantaneous  frequency  sweep. 
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Modelocking 


The  primary  source  for  ultrashort  optical  pulses  is  the  modelcked  laser 
oscillator  [12].  When  many  longitudinal  modes  of  the  laser  cavity  (-186,000  in 
this  case)  are  modelocked  and  equally  spaced  in  frequency  they  combine  to 
generate  short  pulses.  Another  way  to  picture  modelocking  is  a  single  pulse 
propagating  around  the  laser  cavity  must  have  the  same  shape  after  one  round 
trip;  this  is  the  definition  of  a  soliton.  Further  discussion  of  the  soliton  solution 
may  be  found  in  Appendix  A.  Derivation  of  the  Non-Linear  Schrodinger 
Equation.  The  output  of  a  modelocked  laser  is  a  train  of  pulses  at  an  interval 
equal  to  the  round-trip  time  of  the  cavity.  This  concept  was  introduced  earlier 
with  Equation  (1.2)  and  (1.3).  The  shortest  possible  pulse  duration  is  limited  by 
the  inverse  gain  bandwidth  [7]. 

There  are  active  and  passive  methods  to  achieve  modelocked  lasing. 

Active  modelocking  involves  modulating  the  intracavity  loss  or  phase,  while 
passive  modelocking  involves  using  intracavity,  intensity  dependent,  dispersion 
control  elements  that  favor  pulsed  over  CW  operation.  They  cause  cavity  modes 
that  are  above  threshold  to  have  the  same  phase,  and  lock  together  [12].  This 
section  will  focus  on  the  two  passive  methods  for  modelocking  the  Cr:LiSAF 
laser:  the  Kerr-lens  effect  and  a  saturable  absorber  mirror  (SAM). 

The  Kerr  effect  occurs  in  non-linear  crystals,  where  the  index  of  refraction 
depends  on  beam  intensity.  Short  pulses  with  high  peak  intensities  are  well 
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suited  to  induce  the  Kerr  effect.  As  the  intensity  increases,  so  does  the  index; 


the  crystal  acts  like  a  lens,  focusing  the  beam  (Figure  8). 


When  attempting  to  modelock  a  laser  it  is  possible  to  take  advantage  of 
the  Kerr-lens  effect  by  inserting  an  aperture  in  the  beam.  If  the  laser  is 
modelocked,  the  pulses  have  high  peak  intensity,  and  the  Kerr  effect  causes  the 
crystal  to  focus  the  beam.  As  the  beam  is  focused,  the  waist  shifts  from  its  CW 
location  to  a  different  modelocked  location.  If  an  aperture  is  placed  at  the 
modelocked  waist  location,  a  pulse  will  experience  less  loss  than  a  CW  beam. 

By  passing  the  pulses,  the  lens  effect  increases  in  the  crystal  and  the  pulse 
experiences  more  gain  in  each  round-trip,  driving  the  laser  into  pulsed 
modelocked  operation. 

The  other  passive  modelocking  element  used  in  the  Cr:LiSAF  cavity  is  a 
SAM.  A  simple  diagram  is  shown  in  Figure  9.  Each  SAM  must  be  carefully 
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designed  for  a  specific  application.  At  the  far  side,  there  is  a  Bragg-mirror  on  a 


semiconductor  wafer  substrate.  This  is  covered  with  an  absorber  layer  and  then 
a  sophisticated  top  film  system,  which  determines  saturable  loss  [17].  At  high 
intensities,  the  absorber  is  saturated  and  becomes  transparent  to  pulses,  allowing 
them  to  pass  through  to  the  Bragg-mirror,  where  they  are  reflected  back  into  the 
cavity.  At  low  intensities,  the  absorber  is  not  saturated,  and  removes  more  of 
the  incident  energy  from  the  cavity.  This  SAM  loss  mechanism  results  in  phase 
locking  that  causes  pulses  in  the  cavity. 

saturable  absorber  substrate 

<4 

100%  reflector 

Figure  9.  Saturable  Absorber  Mirror  [17]. 

Modelocking  with  the  Kerr  effect  and  SAM  can  only  be  achieved  if  the 
positive  GVD  from  the  laser  crystal  and  the  positive  chirp  from  SPM  are 
corrected  with  negative  GDD  mirrors.  White-light  interferometry  was  used  to 
measure  GDD  in  these  laser  cavity  mirrors. 
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White  Light  Interferometry 


Cr:LiSAF  introduces  frequency-dependent  positive  dispersion  into  the 
laser  cavity.  This  GDD  is  partially  a  result  of  the  ultrafast  SPM  caused  by  the 
optical  Kerr  effect.  Chirped  dielectric  mirrors  with  negative  GDD  balance  the 
Kerr  effect  in  broad  bandwidth  pulses.  A  standard  Michelson  interferometer  and 
a  white-light  source  can  be  used  to  measure  the  dispersive  properties  of  both 
transmissive  and  reflective  optical  components.  The  goal  is  to  measure  the  GDD 
of  the  mirrors  used  in  the  Cr:LiSAF  laser  cavity  in  order  to  determine  if  the  net 
GVD  is  negative.  This  will  provide  the  information  needed  to  choose  mirrors 
with  the  right  amount  of  negative  GDD  in  order  to  modelock  the  laser  and 
generate  terahertz  radiation. 

White  Light 

White  light  is  the  best  choice  to  measure  laser  cavity  mirrors  because  it 
mimics  the  broad  bandwidth  and  short  coherence  length  of  femtosecond  pulses. 

A  Michelson  interferometer  uses  a  beam  splitter  and  mirrors  to  create  an 
interference  fringe  pattern.  Fringes  occur  when  two  or  more  light  waves  yield  a 
resulting  irradiance  that,  deviates  from  the  sum  of  the  component  irradiances 
[18].  If  two  beams  interfere  and  create  fringes,  they  consist  of  nearly  the  same 
frequency.  However,  if  the  beams  are  both  white  light.,  which  contains  many 
frequencies,  then  the  reds  will  interfere  with  only  the  reds,  the  blues  with  only 
the  blues,  etc.,  producing  fringes  as  well. 
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As  a  quick  review,  here  are  some  characteristics  of  white  light,  and  how  it 
lends  itself  to  this  application.  White  light  is  essentially  incoherent,  like  sunlight 
—  unless  one  looks  on  a  femtosecond  time-scale.  Its  limited  temporal  coherence 
is  very  similar  to  ultrafast  pulses.  A  white  light  source  emits  many  wavelengths 
(temporal  incoherence),  and  the  phase  of  the  light  emitted  from  different  regions 
of  the  source  are  not  correlated  (spatial  incoherence).  The  coherence  time,  rQ ,  is 


r 


o 


1 

~Ku 


X 


AX 


(2.1) 


where 


\2 

AX  =  — Av .  (2.2) 

c 

The  longer  the  coherence  time,  the  greater  the  temporal  coherence  of  the  source 
[16].  Another  way  to  describe  coherence  time  is  the  interval  over  which  the  light 
wave  resembles  a  sinusoid.  For  white  light,  rQ  &  1CT14s,  which  is  very  short 
compared  to  a  standard  He-Ne  laser,  where  rQ  «  10  ;;s. 

The  coherence  length,  l  ,  is 


l 


cr„ 


X1 

AX 


(2.3) 


Fringes  can  appear  on  the  detector  when  the  difference  between  the  beam  paths 
to  each  mirror  is  less  than  half  the  coherence  length, 

(.>2(4-0-  (2-4) 

Here,  l  and  Z9  are  the  distances  from  the  beam  splitter  to  each  mirror, 
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respectively.  Additionally,  the  total  phase  difference  must  be  an  integer  multiple 
of  7T ,  and  less  than  the  coherence  length  (2.4).  For  white  light,  l  «  2xlCT4cm, 
which  is  much  shorter  than  a  He-Ne  laser,  where  l  ps  3x10'  cm. 

C 

Interference  Theory 

In  this  experiment,  a  white  light  source  is  directed  into  a  Michelson 
interferometer  (Figure  10).  When  it  reaches  the  beamsplitter,  half  of  the  light 
passes  straight  through  to  the  first  mirror  and  half  of  the  light  is  reflected  90°  to 
the  second  mirror.  After  reflection  on  the  mirrors,  the  beamsplitter  directs  the 
combined  beams  towards  a  detector.  Where  the  reflected  beams  overlap, 
interference  fringes  are  expected  when  the  path  lengths  are  the  same  and  they 
meet  the  criteria  described  in  the  previous  section. 

The  fringe  visibility,  V,  can  be  described  quantitatively  as 


V 


I  -I 

max _ mm 

I  + 1  . 

max  mm 


(2.5) 


The  intensities  in  (2.5)  correspond  to  a  maximum  and  adjacent  minimum  in  the 
fringe  system.  The  greatest  contrast  between  the  dark  and  light  fringes  is  when 
V=1  [18].  Since  the  visibility  depends  on  the  delay  between  the  beams,  r,  it 
can  also  be  described  to  the  first  order  as 


when  t  <  tq  ,  and  zero  otherwise 


v(t)=1-~,  (2.6) 

T« 

[16].  Coherence  time,  r(. ,  is  defined  in  (2.1). 
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Spectrometer 

O  CCD 


Figure  10.  Michelson  interferometer  (MI) 


Arm  2 


Arm  1 


Figure  11.  Simple  Michelson  interferometer 
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Figure  11  shows  a  simple  diagram  of  a  Michelson  interferometer.  A 
sample  mirror  is  placed  in  Arm  2,  and  a  known  dispersive  window  is  placed  in 
Arm  1  with  the  reference  mirror.  This  section  explains  the  theory  behind 
finding  the  GDD  of  the  mirror  in  Arm  2. 

Beginning  with  an  electric  field  E0  incident  on  a  Michelson  interferometer, 
the  electric  field  is  split  into  two  beams.  In  Arm  1,  the  electric  field  after 
traversing  the  arm  and  being  reflected  back  is 

Ei  =  PiE0e~ikw2Le~ik2d ,  (2.7) 

and  the  electric  field  in  Arm  2  after  reflection  is 

E2  =  p2Eoe~‘k2d' .  (2.8) 

The  reflectance  at  each  mirror  is  p{  and  p.2 ;  d  is  the  distance  from  the  beam 
splitter  to  the  first  mirror;  d!  is  the  distance  from  the  beam  splitter  to  the 
second  mirror;  and  L  is  the  thickness  of  the  CaF2  window.  Here,  k  describes 
propagation  through  air,  and  k  describes  propagation  through  the  window. 
Combining  the  field  in  each  arm  yields  the  total  electric  field  at  the  detector 


E  =  E  +  E  . 

T  1  2 


(2.9) 


The  intensity  of  the  combined  beams  is 


I  = 


E„ 


E1  +  E2 


=  E  E  +  E  E  +  E  E  +  E  E 
1  1  1  2  2  1  2  1  1  12 


(2.10) 


Substitute  (2.7)  and  (2.8)  into  (2.10),  and 


/  = 


E, 
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( ~  77!  —ik2d'  \  f  *  77!*  2 L  ik2d  \  ■  (  j~\  — —ik2d  j~i*  ik2d'  \  /r>  i  i  \ 

[P2E  0e  ){PlE0e  e  )+  {PlE0e  e  ){P2E0e  )■  (2J1) 
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To  include  the  phase  change  at  the  of  the  mirror,  let 


(2.12) 


where  p  is  real  because  the  reflectance  of  a  metal  mirror  does  not  cause  a 
frequency  dependant  phase  change,  but  the  reflectance  of  the  dispersive  mirror, 
p  ,  has  a  complex  phase  component.  After  substituting  (2.12)  into  (2.11),  and 
combining  exponentials, 


i2  |  i2  |  1 2  |  i2 


1  =  E0  Pi  +  Eo  P,  I 


I  n  2  -»0„  -ik2(d-d')  ik,.2L 

+  E0  P2  P1  e  Me  K  ’e  w 

I  n  2  ik2(d-d')  —ik,.,2L 

+  EQ  p1  p2  e  Me  y  'e  w 


(2.13) 


Let- A  =  l(d'  —  d^,  the  difference  in  round-trip  path  length,  and  normalize  the 


incident  intensity,  so 


/  2  I  2  I  -;(fcA+<)>  -2/tyL)  i(kA+4> u-2kw  L\ 

— -  =  pl  +  p2  +  p1  p2  e  v  +  e  v  ’ 


(2.14) 


Recall  that 


cosd  = 


i6  .  —id 

e  +  e 


(2.15) 


Now  the  normalized  intensity  becomes 


=  p\  +K|  +2|di||p2|c°s(A'A+^_^2jL)' 


(2.16) 
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If  Arm  1  is  blocked, 


/  oc 


Pi 


E. 


,  and  if  Arm  2  is  blocked, 


I  oc 


P 2 


£ 


.  To 


properly  normalize  the  intensity  of  the  fringes, 


E 


K  2  V  P 1  \  P  2 


(2.17) 


where  the  reference  is  twice  the  square-root  of  the  product  of  each  arm. 

k  =  =27 urn  =  —  n  a  =  —  and  u  =  —  n 

Recall  C  when  .  The  index  of 

refraction,  n,  is  unity  in  air,  but  frequency  dependant  in  general.  Define  the 

argument  of  the  cosine  in  (2.16)  as 

T(u;)=  kA  +  (f)M-kw2L 

The  first  term  of  (2.18)  relates  to  the  field  as  it  travels  through  air, 


(2.18) 


kA  =  —  A. 


(2.19) 


The  third  term  of  (2.18)  relates  to  the  field  as  it  passes  through  the  CaF2 


window, 


*w  =  4 d= 

Using  a  Taylor  series  expansion  of  kw  (u^  at  uQ, 


(2.20) 
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(2.21) 
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where  is  the  group  velocity.  The  second  term  of  (2.18)  relates  to  the  phase 
change  at  the  mirror.  Similar  to  (2.21),  a  Taylor  series  expansion  of  (f>M  at  lo(.  is 


K  (W)  =  K  M- 


d(t>M  (  \  ,  1  d^M  (  \2 

- —  lu-Un)-\ - —  [UO  —  UOn). 

duo  '  oJ  2  duo2  V  oJ 


(2.22) 


(Note  that  (2.22)  is  identical  to  (6.1)  in  Appendix  B).  Combining  equations 
(2.19),  (2.21),  and  (2.22), 


y(uo)=  -A 
v  '  c 

+^V  (%)+ 


d(t>M  (  \  ,  1  d^M  (  f 

den  V  °'  2  duo2  V  °' 


(2.23) 
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c  PC 
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The  term  - —  is  the  GDD  of  the  mirror.  Rearranging  terms  yields 

duo2 


GDD , 


*(")=  — g rDwindow  (w-%) 


den  z/  c 


(cn-cn0) 


(2.24) 
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Recall  that 


(lo  —  uo{)^=  (27 rccr  —  27tc(t0  ) 
=  (27rc)(a-a0) 


(2.25) 


Rewriting  (2.24)  in  the  form  of  a  standard  quadratic  equation,  Ax2  +  Bx  +  C  , 
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/  \  /  \2  GTAD  /  \2 

*W=(H  -^-L  GVD ^  (<7-<7„) 
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da;  z;  c 
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(a~a«) 


+  -A  +  4»G)_2i  n("o) 


The  coefficient  of  the  first  quadratic  term  is 


(2.26) 


/  \2 

T  =  (2ttcJ  - -  L  GVDv 


(2.27) 


Rearranged, 


GDD  = 

Mirror 


-  +  2L  G  VD  , 

\2  Window 


(2.28) 


The  previous  equation  shows  that  the  unknown  group-delay  dispersion  of 
a  mirror  can  be  found  using  this  A  coefficient.  The  theoretical  GVD  of  the  CaF2 
window  may  be  calculated  from  the  Sellmeier  equation  with  coefficients  obtained 
from  SCHOTT  [19]: 


n2(A)=  1  + 


.6188140A2 


.4198937A2 


3.426229A2 


_ I _ | _ (2  29) 

A2  -.002759866  A2  -.01061251  A2  -1068.123’ 


Multiplying  - — -  by  twice  the  thickness  of  the  window  accounts  for  the 


dispersion  that  occurs  as  the  beam  passes  through  the  window  on  its  way  to  the 


mirror  and  again  after  reflection.  It  also  converts  the  window’s  GVD  with 
dimensions  of  time2/length  into  GDD  with  dimensions  of  time2.  At  .846  pm,  256 
fs2  is  the  theoretical  GDD  for  two  passes  through  5  mm  of  CaF2. 
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Curve  Fitting 


The  curve-fit  model  uses  (2.26) 

y  (o ^  =  cos  A  (o  —  oQ  J  —  P  +  S  ,  (2.30) 

where  P  represents  a  phase  factor,  and  S  is  a  vertical  shift.  Each  fit  yields  an  A 
coefficient  which  is  used  in  (2.28)  to  find  GDDMinor  at  the  central  delay 
wavelength  corresponding  to  of) . 

Window  Effect 

Curve-fitting  was  improved  by  using  the  CaF2  window.  The  narrow 
stopband  of  chirped  mirrors  leaves  little  room  for  the  number  of  interference 
fringes  needed  to  apply  the  curve-fit-  model.  A  theoretical  plot  of  (2.30)  using 
(2.27)  with  a  negative  GDD  mirror  alone  in  the  MI  arm  is  shown  in  Figure  12a; 
the  next  plot,  Figure  12b,  shows  the  effect  of  placing  a  dispersive  media  in  the 
reference  arm.  As  the  A  coefficient  become  more  negative,  the  number  of  fringes 
increases.  One  can  see  that  the  second  plot  gives  the  curve-fit  more  information 
to  fit  so  the  coefficients  are  more  accurate. 
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negative  GDD  mirror  with  a  5  mm  CaF2  window  in  the  reference  arm. 


FROG 


In  addition  to  finding  GDD  with  the  Michelson  interferometer,  a  second 
method  was  employed  to  corroborate  the  results.  It  is  a  short-pulse 
measurement  technique  known  as  frequency-resolved  optical-gating  or  FROG.  It 
uses  an  SHG  crystal  and  a  CCD  camera  to  give  near  instantaneous  feedback  on 
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wavelength,  pulse  width,  intensity,  phase  vs.  time,  spatial  profile,  and  spatial 
chirp  [20].  The  Grenouille  FROG,  when  paired  with  the  QuickFrog  software 
package,  measures  ultra-short  laser  pulses  directly,  and  provides  both  spatial  and 
temporal  parameters.  The  Grenouille  measures  the  actual  laser  pulse  width  and 
the  FWHM  value  for  r  .  Although  it  does  measure  spatial  chirp,  one  drawback 
to  the  Grenouille  is  that  is  does  not  provide  any  information  on  the  temporal 
chirp  of  the  pulse.  The  basic  design  of  a  Grenouille  is  shown  in  Figure  13. 
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Figure  13.  Grenouille  FROG  [20]. 


The  Grenouille  FROG  was  used  to  measure  the  change  in  pulse  width  after 
multiple  bounces  between  two  mirrors.  Using  an  equation  similar  the  one  used 
to  find  the  waist  of  a  Gaussian  beam,  the  pulse  train  diverges  like  a  Gaussian 
beam. 
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where  r  is  the  final  pulse  width,  r  is  the  initial  pulse  width,  and  b  is  the 
number  of  bounces  the  beam  makes  on  the  mirrors.  The  dispersion  length 


L  =  — d - .  (2.32) 

9  2  GDD 

Rearranging  (2.31)  gives  just  one  data  point  for  the  GDD  at  the  wavelength  of 
the  TkSapphire  laser  (800  nm) 


GDDU. 

Mirror 


2b 


(2.33) 


Cavity  Dispersion 

Using  the  mirror  measurements  from  the  MI  and  the  FROG,  it  is  trivial 
to  analyze  the  total  cavity  dispersion.  Knox  experimentally  verified  that  the 
dispersion  in  a  cavity  can  simply  be  summed  [21].  Diels  supports  this  result, 
showing  that  linear  optical  elements  may  be  characterized  by  a  complex  optical 
transfer  function 


H(w}=R(wy*(u\  (2.34) 

where  R  is  the  amplitude  response  and  4/  is  the  phase  delay  (2.24) 
experienced  by  a  spectral  component  of  frequency  u .  For  a  sequence  of  rn 
optical  elements,  the  total  transfer  function  is  just  the  product  of  the  individual 
transfer  functions 


H  (u)=  Y\Hj  (u)=  n R.  (-)  exp  (-) 

j= i  vj=1  i=1 


(2.35) 
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The  exponential  term  contains  the  quadratic  phase  response  (GDD),  so  each 
optical  element  in  the  cavity  may  simply  be  summed  to  determine  the  total 
GDD  for  the  system  [10]. 

Now  that  the  groundwork  has  been  laid  for  the  Cr:LiSAF  laser  design, 
and  the  necessity  of  measuring  dispersion  in  the  laser  cavity  mirrors  with  white- 
light  interferometry  has  been  explained,  it’s  time  to  expand  on  the  details  of  the 
experiment . 


30 


III.  Methodology 


The  methodology  section  encompasses  experimental  techniques  used  to 
locate  fringes  with  the  Michelson  interferometer  and  the  details  of  data  collection 
and  analysis. 

Finding  Fringes 

In  order  to  find  fringes,  each  mirror  was  placed  in  the  sample  arm,  while 
the  reference  arm  contained  a  metal  reference  mirror  and  a  CaF2  window.  A 
complex  positioning  process  was  required  to  locate  fringes  within  the  delay  arm’s 
limited  range  of  movement.  The  device’s  original  mirror  had  been  attached  to 
the  front  surface  of  the  holder  with  a  mask.  In  order  to  reproduce  the  location 
of  the  reflective  surface,  while  securing  the  sample  in  the  holder,  the  holder  was 
suspended  above  the  bench  with  two  identical  spacers.  When  the  sample  mirror 
was  placed  into  the  holder  facedown,  it  protruded  the  same  distance  as  the 
original  mirror’s  thickness  (Figure  14).  This  process  meant  that  fringes  could  be 
found  within  the  range  of  the  delay  arm  when  the  sample  was  in  the 
interferometer. 
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Mask 


Sample 


Figure  14.  Sample  holder 


Locating  fringes  followed  the  simple  process  of  aligning  the  reference 
mirror  with  the  sample  mirror,  then  adjusting  the  delay  arm.  To  align  the 
mirrors,  a  crosshair  was  placed  before  the  beam  splitter  and  the  reference  mirror 
was  moved  until  both  crosshair  reflections  overlapped.  The  crosshair  was 
removed  and  the  delay  arm  was  scanned  with  small  movements  until  fringes 
were  located. 

It  was  important  to  make  small  adjustments  when  looking  for  fringes. 

The  fringes  are  often  very  thin,  and  the  range  of  zero-delay  very  short  (±  0.01 
mm).  The  eye  cannot  see  them  with  continuous  movement  of  the  delay  arm 
because  it  integrates  the  intensity  of  the  fringes  during  their  brief  appearance.  It 
is  necessary  to  pause  between  each  small  adjustment,  to  give  the  eye  an 
opportunity  to  see  the  fringes  if  they  exist  at  that  location.  It  is  easy  to  make 
small  adjustments  on  this  interferometer  because  the  delay  arm  knob  is  geared 
with  a  5:1  ratio. 
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Once  fringes  are  located,  minute  tip  and  tilt  adjustments  to  the  reference 


mirror  change  parallel  fringes  into  a  circular  pattern.  In  this  arrangement,  the 
GDD  of  the  sample  mirror  can  be  found  [20].  An  aperture  is  placed  in  the 
output  beam,  and  stopped  down  so  the  spectrometer  sees  only  the  large  center 
fringe. 


Figure  15.  Example  fringes  a.  Linear  [22]  b.  Circular. 


Data  Collection 

When  properly  adjusted  into  a  circular  fringe  pattern,  the  range  of  zero- 
delay  is  no  longer  very  short.  Due  to  the  large  amount  of  GDD  in  these  mirrors, 
fringes  may  be  visible  ±  0.15  mm  on  the  delay  arm.  Scanning  through  this  range 
moves  the  location  of  zero-delay  across  the  frequencies  in  the  white-light.  Figure 
16  shows  zero-delay  at  different  wavelengths  for  three  positions  of  the  delay  arm. 
Most  data  was  collected  at  0.02  mm  intervals  over  the  range  where  fringes  were 
visible.  For  each  sample,  two  reference  files  are  collected.  The  first  is  taken 
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with  the  sample  arm  blocked  and  the  second  is  taken  with  the  reference  arm 


blocked.  The  square  root  of  the  product  of  these  two  arms  is  used  as  the 
reference  to  normalize  each  scan.  This  process  is  described  in  Equation  (2.17). 
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Wavelength  (A) 


Wavelength  (A) 


Wavelength  (A) 

Figure  16.  Example  data  showing  the  zero-delay  wavelength  as  the  delay  arm  is 

scanned 
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Analysis  Process 


It  is  the  regions  of  zero-delay  where  the  theoretical  model  is  fit  to  the 
experimental  data.  The  A-coefficient  from  the  fit  (2.30)  is  used  to  calculate  the 
total  GDD  (2.28)  for  that  wavelength.  This  curve-fitting  process  has  been 
automated  in  MATLAB.  The  program  is  included  in  Appendix  C.  Datacruncher 
m-file.  To  summarize  the  procedure:  the  sample  data  is  normalized  using  the 
reference  data;  the  zero-order  fringe  is  identified;  and  the  model  is  fit  to  a  region 
that  includes  approximately  four  fringes  on  either  side  of  the  delay  (Figure  17). 
The  A-coefficient  and  a  are  saved  to  an  Excel  workbook,  where  er()  is  converted 
to  wavelength  and  A’s  are  used  to  find  the  total  GDD.  Since  all  the 
experimental  data  was  collected  with  the  CaF2  window  in  the  reference  arm,  its 
GDD  must  be  removed  from  the  total  GDD  to  obtain  the  mirror’s  GDD.  This 
process  was  discussed  in  the  section  on  interference  theory. 

The  spectrometer  was  configured  to  maximize  fringe  visibility  without 
saturating  the  CCD  detector  or  losing  the  signal  entirely.  On  the  bench,  the 
detector  and  aperture  were  physically  adjusted  using  near  real-time  feedback 
from  the  program. 
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Figure  17.  Example  curve- fit  to  experimental  data 

FROG 

A  Grenouille  FROG  was  used  to  confirm  the  MI  results.  Two  identical 
mirrors  were  placed  so  that  Ti:Sapphire  laser  output  made  several  bounces 
before  going  into  the  Grenouille  FROG.  Two  different  sets  of  mirrors  were  each 
examined  with  two  to  eight  bounces  (Figure  18).  As  a  reference  (r  ),  the  pulses 
were  measured  with  no  bounces. 


Figure  18.  Six  bounces  before  FROG 
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IV.  Results  and  Analysis 


This  section  presents  the  results  of  measuring  their  GDD  with  the 
Michelson  interferometer  and  the  Grenouille  FROG.  It  addresses  the  offset 
between  these  results,  and  presents  proof  of  modelocked  operation  in  the 
Cr:LiSAF  laser. 

Mirrors 

There  were  six  mirrors  of  interest.  The  Newport  10Q20UF.40  Ultrafast 
Chirped  Negative  GVD  mirrors  are  used  for  multiple  bounces  in  the  cavity.  A 
photograph  of  a  Newport  mirror  is  shown  in  Figure  19.  The  output  coupler, 
made  by  CVI  (PR1-850-99-0525),  is  in  the  position  opposite  the  Newport  mirror. 
This  mirror  is  slightly  less  reflective  than  the  other  mirrors,  allowing  a  portion  of 
the  beam  to  leak  out  of  the  cavity  for  analysis.  The  two  spherical  mirrors  were 
designed  and  manufactured  specifically  for  this  project.  Their  coating  is 
engineered  to  transmit  the  pump  diode  wavelength,  650  nm,  and  reflect  the  846 
nm  emitted  by  the  gain  media.  The  transmission  of  these  mirrors  is  shown  in 
Figure  20.  The  reflectance  (1-T)  is  near  99.8%  over  a  narrow  region  centered  at 
840.  Since  measuring  curved  mirrors  with  the  Michelson  interferometer  is 
experimentally  difficult,  a  third  CVI  mirror  was  used.  This  mirror  has  the  same 
coating  as  the  spherical  mirrors,  but  it’s  flat.  Three  other  mirrors,  not  in  the 
laser  cavity,  are  a  New  Focus  5102  NIR,  Coherent  G0324-005,  and  BATOP 
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SAM  262-IIIa.2.  The  New  Focus  mirrors  were  used  as  a  control  sample,  where 


large  positive  GDD  was  expected.  Like  the  CVI  mirrors,  the  Coherent  case 
contains  two  spherical  mirrors  and  one  flat  mirror  all  with  the  same  coating. 
These  mirrors  are  available  to  focus  the  beam  onto  the  SAM,  if  necessary.  The 
SAM  concept  was  described  in  Chapter  II.  A  photograph  of  a  SAM  is  shown  in 
Figure  19b.  It  is  made  of  a  0.5  in  copper  mount,  with  a  3  mm2  mirror  in  the 
center,  and  it  is  designed  to  have  a  saturable  loss  of  2%  absorption  at  850  nm. 
Even  though  the  SAM  is  so  small  in  diameter  with  respect  to  the  other  mirrors, 
it  was  still  possible  to  measure  fringes  with  the  MI. 


Figure  19.  a.  Newport  mirrors  (1  in)  [23].  b.  BATOP  SAM  (0.5  in)  [18]. 
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Figure  20.  Transmission  of  CVI  TNM2-800-900-0537  99.8%  at  846  nnr. 

MI  Results 

Comparing  the  theoretical  GDD  for  CaF2  with  the  GDD  found  using  this 
experimental  fitting  method  validated  the  analysis  process.  At  the  lasing 
wavelength  of  Cr:LiSAF,  the  theoretical  value  is  256  fs2,  while  the  experimental 
measurement  found  it  to  be  288  fs2.  This  is  within  one  standard  deviation  of  the 
theoretical  value.  Although,  the  difference  is  most  likely  caused  by  the 
beamsplitter  and  compensator  in  the  MI,  since  they  are  each  ~5  mm  thick. 

The  experimental  results  for  all  the  mirrors  are  presented  in  Figure  21. 
The  linear  fit  of  each  data  set  is  shown  along  with  the  R-square  value  indicating 
the  goodness  of  each  fit.  The  first  impression  of  this  plot  is  how  strongly 
dispersive  the  Coherent  and  the  New  Focus  mirrors  are.  The  Coherent  mirror 
has  the  lowest  R-square  value  of  all  the  mirrors.  It  is  likely  that  the  GDD  of 
this  mirror  has  some  modulation  that  does  not  lend  itself  to  a  linear  model. 
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While  the  New  Focus,  Newport,  and  BATOP  mirrors  show  less  modulation  and 


have  better  linear  fits.  Another  noteworthy  plot  is  the  BATOP  SAM.  It  was 
not  possible  to  obtain  good  curve-fits  at  wavelengths  greater  than  700  nm.  The 
equation  of  the  linear  fit  was  used  to  estimate  the  GDD  at  846  nm  where  the 
Cr:LiSAF  lases.  The  same  technique  was  used  for  all  mirrors.  Once  the  value  of 
total  GDD  at  846  nm  was  identified,  the  experimental  value  of  CaF2,  288  fs2, 
was  subtracted  as  described  in  (2.28).  Table  1  shows  the  GDD  measured  for 
each  mirror  at  the  lasing  wavelength. 


Wavelength  (um) 


Figure  21.  Data  showing  the  R2  values  for  linear  fits 
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Table  1.  GDD  results  from  MI 


Mirror 

GDD  at  846  nm  (fs2) 

Newport  10Q20UF.40  -  Ultrafast  Chirped  Neg  GVD 

60.7 

New  Focus  5102-NIR 

420.7 

Coherent  G0324-005 

307.7 

CVI  -  TNM2-800-900-0537-  Flat 

-  6.6 

CVI  -  PR1-850-99-0525  -  Output  Coupler 

-  25.1 

BATOP  SAM  262-IIIa.2 

22.54 

FROG  Results 

The  first  two  mirrors  in  Table  1  were  also  measured  with  the  Grenouille 
FROG,  and  these  results  are  shown  in  Table  2.  The  results  from  the  FROG  are 
closer  to  expected  values.  The  second  column  of  Table  2  shows  the  MI  results 
extrapolated  to  800  nm  from  a  linear  fit  to  each  data  set.  The  third  column 
shows  the  offset  between  the  FROG  results  and  the  MI  results.  The  average  of 
the  differences  is  120  fs2.  This  value  will  be  referred  to  as  the  FROG  offset. 


Table  2.  GDD  results  from  FROG 


GDD  from  FROG 

GDD  from  MI 

Mirror 

at  800  nm  (fs2) 

at  800  nm  (fs2) 

Offset  (fs2) 

Newport  10Q20UF.40 

Ultrafast  Chirped  Neg  GVD 

-77.45 

55.9 

133.35 

New  Focus  5102-NIR 

308.41 

414.27 

105.86 

Table  2  shows  that  the  New  Focus  mirrors  are  much  more  dispersive  than 
the  Newport  mirrors.  It  is  possible  to  see  the  dispersive  effect  as  the  pulse  train 
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mimics  a  diverging  Gaussian  beam  in  Figure  22.  The  experimental  data  is  in 


blue  with  the  theoretical  beam  divergence  in  red.  This  data  was  taken  when  the 
Ti:Sapphire  laser  was  not  very  stable.  Since  the  calculated  GDD  in  (2.31)  is 
very  sensitive  to  the  initial  pulse  width,  r. ,  the  reference  pulse  was  adjusted 
such  that  the  standard  deviation  of  the  results  was  minimized,  within  the  normal 
variance  of  the  input. 

Figure  23  shows  the  GDD  plots  of  these  mirrors.  A  flat  line  is  expected 
from  Equation  (2.31)  since  GDD  is  a  constant  value  for  each  wavelength. 


Figure  22.  Beam  dispersion  from  bounces  between  two  mirrors  a.  Newport  b. 

New  Focus 


Figure  23.  FROG  results  a.  Newport  b.  New  Focus 
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Error  Analysis 


Using  the  MI,  the  standard  deviation  of  the  CaF2  measurements  was 
25.18  fs2,  while  the  standard  deviation  of  the  curve-fitting  process  to  the 
experimental  data  was  25.45x10 "5  pm"1.  Since  there  were  no  known  values  for 
the  mirrors,  the  CaF2  error  may  be  extended  to  all  the  mirrors  measured  with 
the  ML  Using  the  FROG,  the  standard  deviation  is  3.06  fs2for  the  Newport 
mirror,  and  10.5  fs2  for  the  New  Focus  mirror.  This  shows  that  the  MI 
measurements  were  the  largest  source  of  systematic  error. 

After  adding  the  CaF2  window  to  the  beam  with  multiple  bounces  on  the 
Newport  mirrors,  and  observing  its  effect  on  pulsewidth  it  is  certain  that  the 
Newport  mirrors  have  negative  GDD.  So  it  makes  sense  to  consider  shifting  the 
MI  results  by  the  FROG  offset,  120  fs2.  The  cause  of  this  offset  is  unexplained, 
and  further  analysis  is  necessary  to  understand  what  would  cause  such  a 
significant  discrepancy  between  the  MI  and  the  FROG  measurements. 


Table  3.  MI  results  adjusted  for  FROG  offset. 


Mirror 

GDD  at  846  nm  (fs2) 

Newport  10Q20UF.40  -  Ultrafast  Chirped  Neg  GVD 

-59.3 

New  Focus  5102-NIR 

300.7 

Coherent  G0324-005 

187.7 

CVI  -  TNM2-800-900-0537-  Flat 

-145.1 

CVI  -  PR1-850-99-0525  -  Output  Coupler 

-  145.1 

BATOP  SAM  262-IIIa.2 

-  97.46 
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Cavity  Dispersion  Analysis 


The  dispersion  of  each  intracavity  element  has  been  measured  in  isolation, 
and  the  results  have  been  adjusted  using  the  FROG  offset  (Table  3).  The  linear 
cavity  in  Figure  2  uses  the  output  coupler  and  four  CVI  mirrors,  two  Newport 
mirrors,  and  a  SAM.  One  may  simply  combine  their  GDD  to  find  the  dispersion 
of  the  entire  resonator.  Using  (2.35),  the  sum  of  the  crystal  and  the  dispersive 
elements  in  the  cavity  is  -1544.62  fs2. 

—97.46  +  7  x  -145.1  +  12  x  -59-3  +  2  x  75  =  -1544.62  (4.1) 

SAM  CVI  Newport  Cr:LiSAF 

It  was  observed  that  twelve  bounces  between  the  Newport  mirrors  were 
necessary  to  generate  enough  negative  GDD  to  saturate  the  SAM  and  achieve 
CW  modelocking.  Eight  bounces,  resulting  in  -1307.54  fs2  was  not  enough 
negative  GDD  for  this  SAM.  BATOP  no  longer  makes  this  product  because 
they  determined  that  it  is  not  an  effective  device.  So  while  it  is  possible  to 
modelock  the  Cr:LiSAF  laser  using  this  saturable  absorber  mirror,  the  multiple 
bounces,  and  large  amount  of  negative  GDD  may  not  be  necessary  with  a  more 
effective  SAM.  It  is  also  important  to  note  that  if  the  results  from  the 
Michelson  interferometer  (Table  1)  are  not  adjusted  with  the  FROG  offset,  the 
GDD  of  the  modelocked  cavity  is  855.38  fs2.  This  is  far  from  the  negative  value 
that  we  know  is  required  for  modelocking,  and  it  supports  the  argument  for 
using  the  FROG  offset. 
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Modelocked  Results 


Modelocking  was  achieved  with  the  cavity  design  described  in  the 
previous  section  (Figure  2).  The  overall  CW  profile  is  shown  in  Figure  24. 
Figure  25  zooms  in  to  show  a  narrow  region  of  the  modes  in  the  time  domain. 
Figure  26  shows  the  Fourier  transform  into  the  frequency  domain,  where  the 
mode  spacing  is  84.8  MHz,  from  v  =  c/2T .  This  success  validates  using  the 
FROG  offset,  because  without  it,  the  cavity  dispersion  measured  with  the  MI 
would  be  too  positive  to  achieve  modelocking. 
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Figure  24.  Q-switched  modelocking  in  the  time  domain. 


Figure  25.  Zoom  of  time  domain  showing  pulses. 
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Figure  26.  Frequency  domain  showing  modes. 
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V.  Summary  and  Suggestions 


This  chapter  discusses  the  conclusions  that  may  be  drawn  from  the 
results  and  explains  the  impact  of  this  research.  It  also  includes  suggestions  for 
further  research  on  this  project. 

Conclusions  of  Research 

This  research  met  the  goal  of  measuring  the  dispersion  in  laser  cavity 
mirrors  using  white  light  interferometry.  That  task  was  the  bulk  of  the  effort, 
but  near  the  end,  a  cavity  built  with  the  measured  mirrors  was  successfully 
modelocked.  This  validates  the  work  and  shows  that  this  project  is  ready  to 
move  into  the  final  phase  of  generating  terahertz  radiation. 

Impact  of  Research 

The  literature  review  turned  up  papers  that  used  interferometry  to 
measure  dispersive  media  like  glass  [24],  They  mention  that  the  theory  also 
extends  to  reflective  media,  but  none  demonstrate  it  experimentally.  So, 
measuring  dispersive  mirrors  with  this  experiment  was  noteworthy. 

The  most  significant  impact  of  this  research  was  the  development  of  a 
new  technique  to  measure  GDD  in  mirrors  with  a  narrow  stop  band.  After  an 
extensive  literature  search,  no  published  work  was  found  that  uses  a  dispersive 
window  to  increase  the  number  of  spectral  fringes  for  chirped  mirrors. 
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Another  important  contribution  was  a  careful  derivation  of  the 


mathematical  theory  supporting  this  technique.  Most  papers  gloss  over  the 
details,  so  this  was  not  a  trivial  undertaking.  It  also  resulted  in  identifying  an 
error  in  a  published  paper  [24],  This  discovery  is  expanded  on  in  Appendix  B. 
Derivation  of  GDD  Equation. 

Suggestions  for  Further  Research 

The  FROG  offset  needs  further  study.  More  mirrors  should  be  analyzed 
using  the  FROG  and  MI  techniques  to  see  if  the  offset  is  consistent.  The  MI 
setup  should  be  systematically  modified  to  identify  the  source  of  this  offset. 
Different  dispersive  windows,  mirrors  with  wide  stop-bands,  alternate  white-light 
sources  are  all  possible  changes.  Overall,  the  technique  is  effective;  it  is  just  a 
matter  of  finding  the  right  components  to  get  measurements  that  are  more 
accurate.  Fortunately,  the  laser  is  not  very  sensitive  to  the  magnitude  of 
negative  dispersion  in  the  cavity,  so  modelocked  operation  was  still  achieved. 

Measuring  the  cavity  mirrors  and  modelocking  the  laser  completes  the 
first  phase  of  the  larger  project.  The  second  phase  includes  generating  terahertz 
radiation  and  running  on  AA  batteries.  A  photoconduct-ive  switch  may  be  used 
to  generate  terahertz  with  the  modelocked  laser.  Once  the  fundamental  goal  of 
generating  terahertz  has  been  accomplished,  A  A  batteries  can  be  used  to  power 
the  pump  diodes.  Then  the  Cr:LiSAF  laser  will  truly  be  a  portable,  affordable 
source  of  terahertz  radiation. 


49 


Appendix  A.  Derivation  of  the  Nonlinear  Schrodinger  Equation 


This  derivation  follows  the  one  outlined  in  Boyd’s  Nonlinear  Optics  [25]. 


Consider  an  optical  pulse  traveling  in  the  z-direct.ion, 

E  (z,t^=  A  (z,t^e‘^°~  +  c.c. 

Propagate  this  pulse  through  a  dispersive,  nonlinear  medium.  Where 
K  =  nim  (w  )CJ0  /  c  ■  The  pulse  must  obey  the  wave  equation: 

d'2E(z,t)  i  d2D(z,t^j 


(6.1) 


dz 2  c2  dt2 

where  D  contains  both  the  linear  and  nonlinear  contributions  of  the  medium. 
The  Fourier  transforms  of  E  (zj  'j  and  D  (z.  t  j  are: 


(6.2) 


00 

E  (z,t^  =  —  J  E^y-^du 

27T 


oo 

D(z,t)=—J  D{z,u)e~iutdu . 

2,7V 

— oo 

Remember  that  E  (z.  ic  j  and  d(z,uj^  are  related  by: 

D  (z,  =  e  (uj^jE  (z, 

Substitution  of  (6.3)  and  (6.4)  into  the  wave  equation  (6.2),  yields: 

82e(z,uj ) 


(6.3) 


(6.4) 


(6.5) 


8z 2 


+  £  E  (z,uj^=  0  . 


(6.6) 


Focusing  on  the  evolution  of  the  pulse  envelope,  take  the  Fourier  transform  of 


(6.1)  and  substitute  the  result  into  (6.6): 
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oo 

E(z,u)=  f  A  (z, 

— OO 

oo 

=  eik°z  J  A(z,ty^dt 

—oo 

=  a  (z,  to — oj0  yk°z 

where  A  (z, ui j  is  the  Fourier  transform  of  A(z,t^j- 
Now,  the  wave  equation  in  terms  of  A  (z,u  —  is: 

—  Ae^+eL)—  Ae^  =  0 
dz 2  V  V 

2  ikAe^  +  —  -  k2A^z  +  e  L)—Aeik°z  =  0 

0  dz2  0  Wc2 

/  \  d2A 

Since  A  (z,  u>  —  u  )  is  a  slowly  varying  envelope,  drop  the  — -  term.  From 

v  ;  dz 2 

(J1 

electromagnetic  theory,  recall  k2  ^ j  =  e  (u j— .  After  canceling  common 


exponential  terms, 

8A  (z,  u>  —  u>n )  /  \  /  N 

- A_ - +  (k2  -  k2  )A  (z,  u  -  u>Q )  =  0 .  (6.9) 

Because  self-phase  modulation  and  dispersion  are  considered  small  perturbations, 
make  the  approximation  that  k 2  —  k2  =  2 kQ  (k  —  j.  So, 


dA[z,u-  u>)  ,  x  ,  , 

- - - -i  +  2  (k  —  kn)A  (z,  uj  —  un  )  =  0 

dz  V  oy  V’  0/ 

Next,  Taylor  expand  k(u )  about  cjfl  and  include  the  first  nonlinear  term: 


(6.10) 


k(u}  =  kQ  +  n2I  —  +  k' (u  -  )+  -k"(w  -  u;Q  j  .  (6.11) 


51 


Where  I  is  the  intensity 


nc 
~  2t r 


-i  (z,  t )  , 


(6.12) 


k'  is  related  to  the  group  velocity 


da; 


and  k"  is  the  group  velocity  dispersion  (GVD), 


1 

v 


(6.13) 


/."  = 


d2/c  (a;) 


do;2 


1  dv 


v  2  da; 


(6.14) 


Finally,  make  a  change  of  variables  A  =  (u  —  substitute  (6.11)  into  (6.10) 
and  take  the  inverse  Fourier  transform,  A 


A 


dA  (z,  A )  a j  /  \  /  \  k"  „  /  \ 

— ^ - -  —  inj  —  A\z,  Aj—  ik'AA  (z,Aj—i  —  A2A  (z,  Aj  =  0 


(6.15) 


A 


From  the  Fourier  transform  tables: 

{4  (z,  A)}=  A  (z, 

dA(z,t) 

dt 

d2A  (Jm) 


A 


{A.4(z,A)}=! 


(6.16) 


A 


f24'A)}= 


dt2 


Substitution  of  these  Fourier  transforms  yield: 


8 A  (z.  t'j 


dz 


(jJ  ~  /  \ 

inj  —  A{z,tj+  k 


,aA(z’t)l  0. 


dt 


dt 2 


(6.17) 
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Use  a  Galilean  Boost  transformation,  r  =  t - ,  to  change  to  the  reference 

v 


9 


frame  of  the  pulse, 

9  A  (z,  t )  k"  d2A(z,r^)  _  nn2uQ 


+  i- 


=  i  ■ 


A(fir\  A(z,t).  (6-18) 


8z  2  dr 2  27 r 

Equation  (6.18)  is  known  as  the  nonlinear  Schrodinger  equation.  It 

describes  propagation  of  optical  pulses  through  dispersive  media.  The  solutions 


are  solit.ons, 


A  (z,  t)=  Aq  sech 
—2nk" 


/  \  .  k" 

i - z 


T 

T 


0  J 


2  t- 


A 


nn^oT* 


(6.19) 


t  =  pulsewidth 


Note  that  for  AQ  to  be  real  and  positive,  k"  and  n9  must  be  of 
opposite  sign.  Since  n  is  positive  for  Cr:LiSAF,  k"  (GDD)  must  be 


negative  in  the  cavity  to  maintain  pulsed  lasing.  The  required 


magnitude  of  this  GDD  must  be  determined  experimentally. 
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Appendix  B.  Derivation  of  GDD  Equation 


When  curved  fringes  are  observed,  the  curvature  is  quadratic  [22],  The 
quadratic  phase  displacement  may  be  written  as: 


A (j)  (uj^=  L  —  (u  —  0J0 )  k'1 


(8.1) 


)! 

Where  L  is  the  length  of  the  medium,  and  kn  is  the  group  velocity  dispersion  of 
the  medium  at  u>  .  This  equation  comes  from  a  Taylor  series  expansion,  where 
the  first  term  is  linear  dispersion,  and  the  second  term  is  group-velocity 


dispersion.  If  the  wavenumber,  a  —  ^  ,  and 


he 

u>  =  2ttv  —  2i r  —  =  2ixhca . 
A 


(8.2) 


Then  substitute  (8.2)  into  (8.1) 


Acj)  (u> )  =  Lk" 2^  7r2  )i\c2  (a  —  aQ  j 

=  2Ln2c2k"  (a  —  a()  ) 

In  this  situation,  h  is  very  small,  and  can  be  overlooked.  Expanding  (8.3), 

A cj)  (u;j=  Aa2  —  Ba  —  C 

=  2Lir2c2kl,a2  -  ALirV^'a^  -  2Lir2c2k"a2 


(8.3) 


(8.4) 


The  coefficients  are  defined  as 
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(8.5) 


A  =  2Lir2c2k" 

B  =  -AL%2c2k"o[} 
C  =  —2L'K2c2k"  a2 


Now,  it’s  possible  to  find  the  central  wavenumber,  a  ,  in  terms  of  the 


coefficients  by  noticing  that  B  contains  2A: 


B  =  -2Aa 


a°  2  A  . 


It  is  known  that  GVD  (or  k")  is  proportional  to  the  dispersion  in  a  material  [11] 

k”  =  ^—  (8.7 

2t rc2  d\2 


The  original  goal  was  to  find  the  dispersion  in  unknown  materials  using  the 
coefficients  from  a  curve  fitted  to  experimental  data.  So  let’s  try  to  manipulate 
these  equations  into  a  useful  form. 

—  =  AX  kn  =  2nc2k"a03  =  A  a  3  =  X  =  (g  < 

d\2  Ao3  0  7 tL  0  ttLsA b  8nLA2 

Now,  to  find  GDD,  just  multiply  the  GVD  by  the  length  of  material,  L: 


GDD  =  k"*  L  =  — - 


V  A ’lL  V  -Z3  \  -  ~(~2Aro)  -■  j 

27TC2  dA  2  27TC2  8-k\A2  16n2c2A2a  3  2  (Vet 


The  GDD  of  a  mirror  may  be  found  using  just  the  A  coefficient 
of  a  quadratic  equation  fit  to  the  experimental  data.  The  central 
wavelength  for  that  GDD  comes  from  (8.6)  using  both  the  A  and  B  coefficients. 


This  equation  differs  from  the  reference  [24]  by  a  factor  of  two.  When  the 
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experimental  result  was  exactly  half  the  theoretical  result  of  a  known  material,  I 
went  to  the  source  of  the  equation  [22]  and  re-derived  (8.4). 
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Appendix  C.  Datacruncher  m-file 


%% 

function  [dat-acomp]  =  ABDat-acruncher3 

datadir  =  '/Users/Al/Desktop/Fringe  Data/20  Dec/CVI  Thick'; 

%input  is  the  directory  where  the  data  is  located 

%all  files  with  extension  *.DAT  will  be  processed 

%outputs  include  flist ,  the  list  if  filenames  in  order  of  processing 

%These  variables  are  declared  here,  so  they  are  global  for  this  function. 

xdat-a  =0; 

ydata  =0; 

xpos  =  0; 

nfiles  =  1; 

function  button_Callback(source,eventdata) 

%This  function  is  called  when  the  continue  button  is  called 
%Rescale  data  set  to  match  selected  'zoomed  in'  area 
xpos  =  xlim; 
close  (fig); 

xposl  =  floor(xpos(l)); 

xpos2  =  floor  (xpos  (2)); 

xdat-a  =  10000. /xdat-a(xposl:xpos2); 

ydata  =  smooth(ydata(xposl:xpos2),'rlowess',.01); 

%This  section  finds  a  guess  at  the  central  wavenumber  by  looking 

%for  the  maximum  of  the  fourier  domain  versus  the  new  variable  x2. 

%x2  is  a  new  variable  that  goes  as  x~2 

%after  scaling  the  x2fft,  the  new  y's  are  resampled 

x2  =  (xdat-a-centerpos).~2; 

x2max  =  max(x2); 

x2min  =  min(x2); 

x2new  =  x2min+(0:(lengt-h(x2)-l))*(x2max-x2min)/(lengt-h(x2)-l); 
x2fft  =  (0:((length(x2)/2)-l))/(x2max-x2min); 
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y2  =  smooth(interpl(x2,ydata,x2new),'rlowess',.l); 

y2fft  =  abs(fft(y2)); 

y2fft  =  y2fft(l:lengt-h(x2fft)); 

%plot(x2fft,y2fft) ; 

[c,iy2max]=  max(y2fft(3:length(x2fft))); 

fguess  =  2*pi*x2fft(iy2max+2); 

end 

function  ok  _  C  allback  ( sour  ce , event  data) 

%This  function  is  called  when  the  OK  butten  is  pushed. 

%The  structure  datacomp  is  saved  to  the  current  directory  as  a  workplace  file. 

save  datacomp. mat  datacomp; 

nfiles  =  nfiles+1; 

close  (fig); 

end 

function  r  edo  _  C  allback  ( sour  ce ,  event  data) 

%This  function  is  called  when  the  redo  button  is  pushed.  It  simply 
%starts  the  data  analysis  over  with  the  current  file, 
close  (fig); 
end 

function  txt  =  myupdatefcn(empt,event_obj) 

%funct-ion  called  when  the  plot  is  clicked 

pos  =  get  (event  _obj,  'Position'); 

centerpos  =  10000/xdata(pos(l)); 

txt  =  {['Wavelength:  ',num2str(xdata(pos(l)))]}; 

line([pos(l)  pos(l)],[0  3]); 

end 

%This  is  the  beginning  of  the  function 

%It  reads  the  list  of  data  files  and  the  reference  file. 

%Then  the  reference  file  is  smoothed.  Before  entering  the  major  while  loop, 
datalist  =  dir(fullfile(datadir,'*.DAT')); 
reflist  =  dir(fullfile(datadir,'*.ref)); 
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refinput  =  importdata(fullfile(datadir,reflist(  l).name)); 
refinput  =  refinput. data; 
refdat-a  =  smooth(refinput(:,2),.01,'rlowess'); 
while(nfiles  <=  length (datalist)) 

%Read  in  the  data  file,  normalize  it  with  the  ref  and  smooth  data. 

datainput  =  importdata(fullfile(datadir,datalist(nfiles)  .name)) ; 

datainput  =  datainput.  data; 

xdata  =  datainput  (:,1); 

ydat-a  =  dat-ainput(:, 2). /refdat-a; 

ydata=smooth(ydata,.01,'rlowess'); 

%Generates  data  plot  and  turns  on  the  cursormode 

fig  =  figure  ('Toolbar', 'figure'); 

plot(ydata); 

ylim([0  4]); 

zoom  on; 

dcm  obj  =  datacursor  mode  (fig); 

set  (dcm_obj  ,'UpdateF  cn',@myupdatefcn) ; 

%  Help  directions  printed  to  screen 

t-ext(0,2.9,'Zoom  to  select  range.  Cursor  to  select  center', 'FontSize', 12); 
t-ext-(0,2.6,'Then  press  continue', 'Font-Size', 12); 

%Generate  continue  push  button  and  wait-  for  user  t-o  to  push  continue 

pbh  =  uicontrol(fig, 'Style', 'pushbutton', 'String', 'Continue', 'Position', [10  10  60 

40], 'Callback', {@button_Callback}); 

uicontrol(pbh); 

uiwait(fig); 

%call  nonlinear  fit  using  the  function  fringe,  fguess  and  the  selected 
%center  point  are  input  into  the  guess. 

[bet-ahat-,resid,J]  =  nlinfit  (xdata, ydata, 'fringe',  [fguess, centerpos, 3, 1,1]); 

%The  fitting  routine  could  be  called  a  second,  or  third  time;  however 
%it  did  not-  seem  t-o  help. 

%[betahat,resid,J]  =  nlinfit- (xdata, ydat-a,  'fringe',  bet-ahat-) ; 
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%betahat 

%betahat,resid,J]  =  nlinfit  (xdata, ydata, 'fringe', betahat) ; 

%betahat 

%nlparci  returns  the  95%  confidence  interval  for  the  fitted  parameters 
betaci  =  nlparci  ( betahat  ,resid,J); 

A(nfiles)  =  betahat  (1); 

sigmaA(nfiles)  =  abs(betaci(l,l)-betaci(l,2))/2; 
sigmaO(nfiles)  =  betahat(2); 

sigma0sigma(nfiles)=abs(betaci(2,l)-betaci(2,2))/2; 
yfit  =  fringe  (betahat,  xdata); 

%create  arrays  of  data  and  store  in  structure  datacomp 
flist.  {nfiles}  =  datalist(nfiles).name; 
datalist.  ( nfiles) .  name 

datacomp  ( nfiles)  =  struct  ( 'filename'  ,datalist  ( nfiles) .  name , 

'Acoefficient' ,  A  ( nfiles)  ,'sigmaA'  ,sigmaA  (nfiles)  ,'sigmaO'  ,sigmaO  ( nfiles)  ,'sigmaOsigm 

a',sigmaOsigma(nfiles)) ; 

plot  (xdata,  ydat-a,'*',xdata,yfit,'+'); 

%generate  ok  and  redo  push  buttons 

pbok  =  uicontrol(fig, 'Style', 'pushbutton', 'String', 'Ok', 'Position', [10  10  60 

40] , '  C  allback ',{  @ok  _  C  allback  }) ; 

uicontrol(pbok); 

pbredo  =  uicontrol(fig, 'Style', 'pushbutton', 'String', 'Retry', 'Position', [10  50  60 

40], 'Callback', {@redo_Callback}); 

uicontrol(  pbredo) ; 

xlabel( 'Wavenumber  ( 1  / microns) ') ; 

ylabel( 'Fringe  Amplitude  (normalized)'); 

title('Fitting  Fringe'); 

legend  ('data',  'Fringe  Fit'); 

uiwait(fig); 

end  %while  loop  end 

end  %function  end 
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Appendix  D.  Fringe  m-file 


function  yhat  =  fringe  (beta, x) 

A  =  beta(l); 
b  =  beta(2); 
c  =  bet.a(3); 
d  =  beta(4); 
f  =  beta(5); 

yhat  =  f+d*cos(A*(x-b).~2+c); 
end 
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